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ABSTRACT 

Despite being the prototype of its class, T Pyx is arguably the most unusual and poorly understood 
recurrent nova. Here, we use radio observations with the VLA to trace the evolution of the ejecta 
over the course of the 2011 outburst of T Pyx. The radio emission, which is consistent with thermal 
emission from the nova ejecta, began rising surprisingly late in the outburst, indicating that the bulk 
of the radio-emitting material was not ejected until ^50-100 days after optical discovery. Given that 
the radio-emitting material must have been ionized, and taking into account the resulting constraints 
on the electron temperature and dumpiness, the high peak flux densities of the radio emission require 
a massive ejection of [2 — 33] x 10~ 5 Mq. This ejecta mass is much higher than previous theoretical 
predictions of outbursts in T Pyx, and implies a very high mass accretion rate onto the white dwarf 
between outbursts. 

Subject headings: white dwarfs — novae, cataclysmic variables — stars: individual (T Pyxidis) - 
radio continuum 



1. INTRODUCTION 

Nova outbursts are the most common stellar explo- 
sions in the universe, and are important so urces of both 
energ y and matter injection into the ISM (|Gehrz et al.l 
1998). Novae are understood to be the result of a ther- 
monuclear runaway on the surface of an accreting white 
dwarf (WD). This runaway begins when the critical pres- 
sure and temperature required for CNO cycle burning of 
hydrogen are rea c hed at the base of th e accreted shell 
(|SchatzmaiJll94ft iStarrfield et al.lll972l) . As long as ac- 
cretion continues, all novae are expected to recur on 
times cales from ~ 10 8 years to less than a year (see 
e.g. lYaron et alJ 12005). This recurrence time depends 
primarily on the mass of the white dwarf and the ac- 
cretion rate from the binary companion; more massive 
white dwarfs have smaller radii and therefore higher sur- 
face gravities, meaning that the critical conditions are 
reached for smaller accreted envelopes. The accretion 
rate controls the time required to accumulate the critical 
mass — higher accretion rate systems accrete the trigger 
mass in less time. 

One of the most important predictions of theoretical 
studies of novae is that m ost systems experi ence a net loss 
in mass during outburst (jYaron et al.ll2005l) , implying se- 
vere difficulties in growing WDs to the Chandrasekhar 
mass so that they explode as SNe la. Only white dwarfs 
accreting material at very high rates from their compan- 
ions are predicted to eject less mass than had been ac- 

t nelson® phy sics.umn.edu 

School of Physics and Astronomy, University of Minnesota, 
115 Church St SE, Minneapolis, MN 55455 

2 Department of Physics and Astronomy, Michigan State Uni- 
versity, East Lansing, MI 48824, USA 

3 National Radio Astronomy Observatory, P.O. Box 0, So- 
corro, NM 87801, USA 

4 Columbia Astrophysics Laboratory, Columbia University, 
New York, NY, USA 

5 Center for Space Science and Technology, University of 
Maryland Baltimore County, Baltimore, MD 21250, USA 

6 CRESST and X-ray Astrophysics Laboratory, NASA/GSFC, 
Greenbelt MD 20771 USA 



creted in the quies cence period prior to outburst (see e.g. 
lYaron et al.l 120051 ). Such high accretion rates are also 
expected to result in relatively short intervals between 
outbursts, with massive white dwarfs experiencing out- 
bursts every decade or so — the observed characteristics 
of recurrent novae (i.e. the subclass of novae with more 
than one recorded outburst). Recurrent novae present 
the best opportunity to determine whether — and in what 
circumstances — accreting WDs can be SN la progenitors, 
by measuring which is larger: the mass ejected in the 
nova outburst or the mass acc reted since t h e last out- 
burst. However, the models of lYaron et al.l (|2005l ) find 
that the difference between the accreted and ejected mass 
is seldom greater than a factor of ~2. Therefore to truly 
constrain whether a WD is growing or shrinking in mass, 
we need accurate and precise measurements of both the 
accretion rate and the ejecta mass. Radio observations 
are key in determining the latter. 

T Pyx is the prototype recurrent nova, with seven out- 
bursts known between 1890 and 1966. However, despite 
its canonical status and more than a century of study, 
many questions about this system remain unanswered. 
T Pyx differs from most other recurrent novae in a num- 
ber of important ways. The outburst evolution is much 
slower, taking longer to both rise to ma ximum light an d 
subsequently fade back to quiescence (|Schaeferi 12010). 
T Pyx is the only recurrent nova w ith an orbital period 
below the period gap (P = 1.83 h; ISchaefer et alj|1992t 
IPatterson etal]|1998t lUthas et all 120101 ) . Masslransfer 
in such low period systems is driven by angular momen- 
tum loss due to gravitational radiation, which typically 
leads to accretion r ates of order 10~ 10 M Q yr" 1 or less 
(jKnigge et al.ll201ll ). yet ev idence abounds for a high ac- 
cretion rate in this system (jSelvelli et al.|[2008l ). Finally, 
the binary system is surrounded by a clumpy remnant 
that is not observed in other recurrent s ystems, extend- 
ing up to 10" from the central binar y (Williams! 119821 : 
IShara et alJfl997t ISchaefer et alJl20loh . 

These unusual characteristics have led to questions 
about the evolutionary state of T Pyx. The short or- 
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Figure 1. Optical and radio light curves of the 2011 T Pyx outburst. Upper panel: V band optical light curve from AAVSO data. Lower 
panel: VLA light curves at frequencies from 1.5 to 33 GHz. We take 2011 April 14 as to, the start of the outburst. 



bital period and low mass of the companion star make it 
difficult to maintain the high accretion rate required to 
power a nova outburst every few decades. To solve this 
problem, authors have suggested that a normal nova out- 
burst in the late 1800s led to a period of enhanced mass 
loss induced by irradiation of the co mpanion star by the 
hot, shell burning w hite dwarf (e.g., iKnigge et al.ll2000l ; 
iSchaefer et aL|[2010h . This higher mass loss rate in turn 
fed a higher accretion rate onto the white dwarf, and led 
to the accumulation of the critical mass required for an 
outburst in just a decade or so. There is evidence that 
the average accretion rate between outbursts has been 
declining over the 20th century, and that the recurrence 
time between novae has been growing, suggesting that 
the current state of the T P yx system is unsusta inable 
for extended periods of time ((Schacfcr et al.ll2010h . 

T Pyx experienced its most recent outburst in April 
2011 — the first since 1966, and thus the first eruption 
for which the astronomical community has been able to 
obtain broad multiwavelength coverage including the ra- 
dio regime. The F-band optical light curve of the out- 
burst, created using data from the American Association 
of Variable Star Observers (AAVSO), is shown in the up- 
per panel of Figure 1. The optical evolution of T Pyx 
during outburst is much slower and more complex than 
that observed in most other recurrent novae. The system 
brightened quickly from its quiescent value of V ~ 15.5 
mag to V ~ 8 mag within the first 2 days of outburst. 
This initial brightening was then followed by a plateau 
phase, during which the optical brightness remained at 
^8 mag until around Day 15. T Pyx then began to 
brighten again, although with some fluctuations, reach- 
ing a maximum brightness of V « 6.5 mag around Day 
30. Over the next two months, the optical light curve 
declined relatively gradually, with some brightness fluc- 



tuations, until around Day ^90 at which time the decline 
suddenly steepened. At the time of completion of this 
work (October 2012, roughly 1.5 years after the outburst) 
the V band magnitude as reported by AAVSO observers 
was ^14.5 mag, indicating that the syst em has not yet 
fully returned to its quiescence flux levels (Schacfcr ct al. 

[Iml). 

Radio emission in novae traces the bulk of the ejected 
material in a relatively simple manner — through ther- 
mal free-free emission. As the thermal radio emission 
transitions from optically thick to optically thin and the 
radio photosphere recedes through the ejecta, the radio 
light curve surveys the entire density profile of the ejected 
mass. In this paper we present multifrequency radio light 
curves of the 2011 nova outburst in T Pyx, obtained with 
the newly-upgraded Karl G. Jansky Very Large Array 
(VLA) of the National Radio Astronomy Observatory, 
and use them to determine the mass of the ejected mate- 
rial. In Section 2, we present our observations and data 
analysis. In Section 3, we describe the evolution of the 
radio emission over the course of the outburst. In Section 
4, we estimate the mass of the ejecta and the timing of 
the ejection. Finally, we present our conclusions in Sec- 
tion 6. Throughout this paper, we assume a distance to 
the system of 4.8 ± 0.5 kpc, derived from light echo ob- 
servations with the Hubble Space Telescope (HST) during 
the 2011 outburst (Sokoloski et al. 2012, in preparation). 
This is consistent with the distance of >4.5 kp c inferred 
from optical spectroscopy bv lShore et al.l (|2011| ) . We also 
take 2011 April 14.0 UT (MJD = 556 65.0) as t , the start 
of the outburst (jWaagan et al .1 [20TTh . 

2. OBSERVATIONS AND DATA REDUCTION 

We obtained sensitive radio observations of T Pyx be- 
tween 2011 April 22 and 2012 September 23 with the 
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Figure 2. A selection of VLA radio spectra over the evolution of the outburst. At each epoch, we over-plot the best-fit single or broken 
power law spectrum in red. A single index quoted in the top-left corner indicates that a single power law is sufficient to fit the data. 
Where two indices are quoted, oti describes the lower frequency data while oli is the fit to the higher frequency data; the break frequency 
is denoted as v^. Where low- and high-frequency data are not simultaneous, we fit them separately and report two indices. 



VLA through programs 11A-263, 11B-230, and 12A-459. 
Over the course of the outburst, the VLA was operated 
in all configurations, and data were obtained in the L (1— 
2 GHz), S (2-4 GHz), C (4-8 GHz), X (8-8.8 GHz), Ku 
(12-18 GHz), K (18-26.5 GHz), and Ka (26.5-40 GHz) 
bands, resulting in coverage from 1-37 GHz. Observa- 
tions were acquired with 2 GHz of bandwidth (the maxi- 
mum available), split between two independently tunable 
1-GHz-wide basebands, with the exception of L-band ob- 
servations, which provide 1 GHz of bandwidth between 
1-2 GHz, and X-band observations, which were still sub- 
ject to the old receivers and therefore only cover 800 MHz 
of bandwidth. The details of our observations are given 
in Table 1. 

At the lower frequencies (L, S, and C bands), the source 
J0921-2618 was used as the gain calibrator, while J0900- 
2808 was used for gain calibration at the higher frequen- 
cies (X, Ku, K, and Ka bands). The absolute flux den- 
sity scale and bandpass were calibrated during each run 
with either 3C138, 3C147, or 3C286. Referenced pointing 
scans were used at Ku, K, and Ka bands to ensure accu- 
rate pointing; pointing solutions were obtained on both 
the flux calibrator and gain calibrator, and the pointing 
solution from the gain calibrator was subsequently ap- 



plied to our observations of T Pyx. Fast switching was 
used for high-frequency calibration, with a cycle time of 
^2 minutes. Data reduction was carried out using stan- 
dard routines in AIPS and CASA. Each receiver band 
was edited and calibrated independently. The calibrated 
data were split into their two basebands and imaged, 
thereby providing two frequency points, with the excep- 
tion of X band (where only 800 MHz of bandwidth is 
obtained, centered at 8.4 GHz) and K band (where the 
two basebands were placed side-by-side and averaged to- 
gether). In most cases, when there was sufficient signal, 
a single iteration of phase-only self-calibration was car- 
ried out, with a solution interval of 1-5 minutes, using 
the clean components from the first imaging run. 

One complication affecting our data was that A con- 
figuration (the most extended VLA configuration) oc- 
curred during the summer months (June- August 2011), 
when the VLA site is subject to poor weather. These 
adverse conditions, combined with the low elevation of 
T Pyx, led to severe decorrelation of many of our A con- 
figuration observations, particularly at high frequencies. 
The K and Ka band data, obtained on 2011 August 20 
and September 17, were in particular very badly decor- 
related. Despite this problem, significant flux density 



4 



was clearly detected; a point source model was used in 
self-calibration to recover the flux. 

In each image, the flux density of T Pyx was measured 
by fitting a gaussian to the imaged source with the tasks 
JMFIT in AIPS and gaussf it in CASA. We record the 
integrated flux density of the gaussian; in most cases, 
there was sufficient signal on T Pyx to allow the width 
of the gaussian to vary slightly, but in cases of low signal- 
to-noise ratio, the width of the gaussian was kept fixed 
at the dimensions of the synthesized beam. Errors were 
estimated by the gaussian fitter, and added in quadra- 
ture with estimated calibration errors of 5% at lower fre- 
quencies (<10 GHz) and 10% at higher frequencies (>10 
GHz). All resulting flux densities and uncertainties are 
presented in Table 1. 

3. RADIO EVOLUTION OF THE 2011 OUTBURST OF 
T PYX 

The VLA light curves, presented in the lower panel of 
Fig 1, show the evolution of the radio emission over the 
course of the outburst. The source was not detected at 
either 33 or 6 GHz during the first observation, on Day 
8, with upper limits on the flux density <0.1 mJy. Sub- 
sequently, on Day 17, T Pyx was detected at 33 GHz at 
a significance of 8a, but was not detected at 6 GHz. In 
the next two epochs (Days 31 and 45-48), we observed 
with additional receiver bands, but the source was not 
detected at any frequency. The most constraining flux 
density upper limits at this time were obtained at high 
frequencies, again with F v < 0.1 mJy. On Day 128, ob- 
servations at 5 and 7 GHz finally resulted in detections 
of the source at the ~0.1 mJy level, and in subsequent 
observations T Pyx brightened at all frequencies, reach- 
ing a maximum flux density of 48 mJy at 37 GHz on Day 
156. The radio flux began to fade after this date, first 
at the highest frequencies, and later at lower frequen- 
cies. By Day 538 (the final epoch of observations) the 
flux density had declined at all frequencies to a level of 
0.5-0.9 mJy. 

At each epoch, we modeled the resulting spectra as 
a single power law of the form S v oc v a (where a is 
the spectral index) , and when flux density measurements 
were obtained at more than 5 frequencies, we also fit a 
broken power law. The fitting procedure was carried out 
using mpf it, a cur ve fitting package implemented in IDL 
(|Markwardtll2009T) . For each spectrum, we compared the 
goodness of fit of the single and broken power law models, 
as measured by \ 2 jv (where v is the number of degrees 
of freedom in the model). We selected the model with 
the lower value of reduced x 2 as the best fit, and plotted 
it as a red line in each panel of Figure 2. During the 
early epochs, the broken power law model provides only 
a marginally better fit to the data than the single power 
law. However at later times, where the turnover at higher 
frequencies is quite clear, the broken power law provides 
a large improvement in the goodness of the fit. 

The best fit spectral indices at each epoch indicate an 
evolution from optically thick emission during the ris- 
ing phase, to optically thin emission at late times. We 
present a selection of spectra that best demonstrates this 
evolution in Figure 2. The spectral indices fit to both the 
low- and high-frequency data are initially steeply rising 
(ai =1.75 ± 0.09, a 2 =1.48 ± 0.06 on Day 129), in- 
dicative of optically thick emission. We note that spec- 



tral index implied by the data on Day 17 (a > 1.14) is 
also steep, and suggests that the source was also opti- 
cally thick at the time of first detection. As the outburst 
evolves the spectrum becomes shallower, beginning first 
at high frequencies on Day 194 , and eventually at low 
frequencies by Day 300. The best fit indices at high fre- 
quencies are between a — 0.3 and 0.8, consistent with the 
values expected for a partially optically thick shell with 
a density profile that decreases in the outward direction. 
By Day 300, the high frequency spectrum has turned 
over (a = -0.23 ± 0.07), consistent with emission from 
a completely optically thin shell. By the epochs of our 
last observations (Days 512 and 528), the spectra have 
indices expected for optically thin gas at all frequencies. 

4. EJECTA MASS AND EJECTION TIMING 

The radio light curves, and the evolution of the radio 
spectra from optically thick to thin, are consistent with 
thermal emission from the expanding nova ejecta. Con- 
sider an ionized shell of material, ejected from the surface 
of the white dwarf shortly after the onset of the nova out- 
burst. The shell is a source of radio flux due to thermal 
free-free emission. At early times the shell is compact 
and dense, and is optically thick at all radio frequencies. 
This yields a steeply rising spectrum with a — 2, as ex- 
pected for the Rayleigh- Jeans tail of a blackbody. At 
this stage, the radio photosphere (i.e. the surface where 
optical depth r = 1) at all frequencies is essentially coin- 
cident with the outer surface of the ejecta, and increases 
in surface area as the ejecta expand. This leads to a 
brightening of the radio emission, and accounts for the 
rising, steep spectrum portion of the radio light curve. 

Some time later, the density in the shell drops suffi- 
ciently that the ejecta start to become optically thin. 
This processes begins at high frequencies first, since 
higher radio frequencies require larger emission measures 
to be optically thick compared to lower frequencies. The 
radio photosphere at such a frequency will recede though 
the ejecta, tracing slower moving ejecta as it goes. This 
leads to a gradual turnover in the radio light curve, 
and a flattening of the spectral index (to a « 0.6 for 
an r~ 2 density profile). This process cascades towards 
lower and lower frequencies as time passes. Once the 
photospheres of all observed frequencies have completely 
receded through the entire ejected shell, the emission 
can be described completely by optically-thin thermal 
bremsstrahlung (a = —0.1). After this, the radio light 
curve simply fades at all frequencies as the ejecta expand 
and the emission measure continues to drop. 

We can estimate the mass of the ejecta using simple 
models that account for the observed radio flux evo- 
lution at all frequen cies in a consistent m anner . Fo l- 
lowing the wor k of iSeaauist k, Palimakal (e.g., 11977ft : 
iHiellming et al.l (e.g.. I1979D . we assume that the ejecta 
can be described as a thick shell of warm (T » 10 4 K) 
material emitting thermal bremsstrahlung radiation. In 
this model, the shell is ejected at to, i.e. the start of the 
nova outburst. The density profile of the shell varies with 
radial distance such that p(r) oc r~ 2 or r~ 3 , with an in- 
ner boundary to the ejecta. The thickness of the ejected 
shell may be due to either a velocity spread in a homol- 
ogous explos ion (a Hubble flow) or a prolonged period 
of mass loss (|Seaauist fc Bodd l2008). The model traces 
the evolution of the photosphere at each frequency as 
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Figure 3. Hubble flow model for the radio emission from the nova 
ejecta, assuming mass ejection at the time of optical discovery. We 
assume a distance to T Pyx of 4.8 kpc, and expansion velocity 
of 3,000 km s _1 . The model plotted here has ejected mass of 
2.5 X 10~ 4 Mq, T e = 10, 000 K, and ti min /u max = 0.9. The match 
to the high flux densities observed at peak brightness is relatively, 
but the model fails to reproduce all other features of the radio 
evolution. 



Figure 4. Same as Figure 3, but allowing the time of mass ejection 
to be an additional free parameter. The formal best fit model has 
t ej = 96 days, M ej = 2.2 X 10~ 4 M , and T e = 67,000 K, and 
"min/"max = 0.2, although we note that the uncertainties on the 
ejection date are large. 



densities indicates that the density in the ejecta, and 
hence the mass, must be high. 



the ejecta expand over time, and depends on the ejected 
mass (M e j), the velocity of the outermost ejecta (w max ), 
the electron temperature of the ejecta (T e ; assumed to be 
constant in time and throughout the ejecta), and the dis- 
tance to the nova (d). The turn-over in the light curves 
also depends on the assumed density profile (we take 
p(r) oc r~ 2 here) and on the thickness of the shell, which 
we parameterize here as the ratio of the velocities in the 
innermost and outermost ejecta (v m m/fmax), assuming 
an instantaneous and homologous explosion. We fit this 
model of expanding thermal ejecta to all observed fre- 
quencies simultaneously with the singular exclusion of 
the lower baseband in L band (~1.3 GHz; this frequency 
mostly yields non-detections ), making use of t he mpf it 
package developed for IDL (|Markwardtl 120091 ) to carry 
out least-squares fitting and obtain the best-fit param- 
eters. Measurements that are not significant at the 3ct 
level are excluded from our model fits. 

The large flux densities reached at peak radio bright- 
ness around Day 200 place important constraints on the 
mass of the ejected shell. Assuming a distance to T Pyx 
of 4.8 ± 0.5 kpc (Sokoloski et al. 2012, in prepara- 
tion), a shell temperature T e of 10,000 K, and a max- 
imum velocity of the ejecta of 3,000 km s _1 (the value 
re ported around t he tim e the radio flux begins to rise 
bv lOsborne eTall (|20lTh ). we find that only models with 
M e j « 2.5 x 10~ 4 M Q can reproduce the peak flux den- 
sities. An example of such a model is shown in Figure 
[3J Irrespective of the quality of the fits at early and late 
times, the requirement that the ejecta must be large and 
optically thick to produce the high observed radio flux 



4.1. Modeling a Delayed Ejection 

The rise of the radio flux is not consistent with the 
primary ejection of material at tQ. At early times while 
the ejecta are optically thick, flux densities are expected 
to rise as approximately t 2 as the ejecta expand. If the 
material responsible for the radio emission had started 
expanding at to, we would have detected brightening ra- 
dio emission at early times (see Figure [3]). Our sensitive 
upper limits rule this out. Despite successfully repro- 
ducing the high-frequency peak flux densities, the model 
shown in Figure [3] fails to reproduce all other details of 
the observed light curve — the rise of the radio flux, the 
spectral evolution, and the late decline phase. The ob- 
served radio rise is steeper, and begins later in the out- 
burst. Although there is a significant detection of flux 
on Day 16, the source is faint (0.248 mJy at 33 GHz) 
and is subsequently undetected at the same frequency in 
the following two epochs. This early detection detection 
may be related to some early period of mass loss (see 
Section 5.3), but it is clearly not related to the primary 
mass ejection traced by the bulk of the radio emission. 
This begins much later, starting with the detections at 5 
and 7 GHz on Day 66. 

We obtain much better model fits to the data if we 
do not require mass to be ejected at the time of optical 
discovery, and instead allow for a delayed mass ejection. 
Including the time of ejection (measured relative to the 
date of discovery, 2011 April 14) as an additional free 
parameter, we find that the best fit models are those 
in which the nova shell is ejected between 50 and 100 
days after optical discovery. An example fit, in which 
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we continue to assume v max =3,000 km s _1 and d = 4.8 
kpc, is shown in Figure 2) The best fit ejected mass is 
M C j = 2.2 x 10~ 4 , the ejection time to is 96 days, the 
electron temperature of the ejecta is T e = 67, 000 K, 
and v m in/«max = 0.2. Not only is this model consistent 
with the non-detections at early times, it provides a bet- 
ter fit to data at late times than the model in Figure 
3, reproducing the steep decline observed in the multi- 
frequency light curve. The relatively high electron tem- 
perature found by the model fit is required to account 
for both the high peak flux densities and the observed 
rate of transition to an optically thin spectrum. 

The fit to the data is shown in more detail in Figure 
[5l which plots the light curve at each frequency in sepa- 
rate panels. The model described above is over-plotted 
as a red line in each panel; it is an excellent fit to the 
light curves at high frequency (>13 GHz). However, at 
lower frequencies the model and the observed light curves 
diverge, with the measured flux densities rising and de- 
clining earlier than the best fit model predicts. Although 
the formal best-fit delay time is 96 days, our detection 
of optically thick emission on Day 66 that brightened 
in subsequent observations suggest that the actual de- 
lay time is probably closer to 50 days. The more grad- 
ual rise observed in the data, and the discrepancies be- 
tween low and high frequencies, suggest that our simple 
model fails to include all the important physics during 
the earliest stages of mass ejection (see Section 5.3 for 
more discussion). Discrepancies between the data and 
the model during the light curve decline might be ex- 
plained by non-isothermal ejecta, which likely exhibit a 
radial temperature gradient. Modeling of non-isothermal 
ejecta is outside the scope of this paper, but has been re- 
quired in the past to explain the radio light curve of the 
nova V1974 Cyg (IHiellming|ll9"96l) . 



4.2. Allowed ranges of the ejecta mass and electron 
temperature 

There are significant degeneracies between the four 
fundamental parameters of the model (M e j, v max , T e ; and 
d), as illustrated in Figure |6l Some of these parameters 
are poorly constrained by the data, for example t> ma x, 
the maximum velocity of the ejec ta. Although we take 
the value 3000 km s -1 reported bv lOsborne et alj f|201 lh 
as our fiducial value, a range of velocities have been 
repor ted in T Pyx during outburst (see e.g. iCatchpoTel 
1969). Similarly, although we favor the distance of 4.8 ± 
0.5 kpc found using light echos from the nebular remnant 
observed with HST (Sokoloski et al. 2012), distance esti- 
mates in the literature range from 2.5 kpc to >5 kpc. In 
order to understand the range of allowed ejected masses 
associated with these uncertainties, we explored a plau- 
sible parameter space of ejecta velocity and distance; for 
each combination we obtained the best-fit ejected mass 
and electron temperature, fitting to all frequencies simul- 
taneously. The results are shown in Figure [6] 

A change in the distance simply moves the model 
light curves up and down in brightness, but it influences 
the derived M e j, because more massive ejecta will have 
higher peak flux densities. Therefore, if T Pyx is located 
at a smaller distance, we find a smaller ejecta mass from 
our radio light curves (top panel of Figure IB])- However, a 
smaller ejecta mass will also become optically thin more 
quickly because the density in the shell is lower. The 
electron temperature covaries with ejecta mass because 
higher temperatures make the ejecta brighter while they 
are optically thick, but also lead to them becoming op- 
tically thin faster. Therefore, M e j and T e are somewhat 
inversely related in then model fits, as seen by comparing 
the two panels in Figure |6l Finally, the expansion veloc- 
ity essentially acts as a "stretch" parameter in the time 
dimension, with slower values of w max producing slower 
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Figure 6. Illustration of the degeneracies in our model fit to the radio light curves of T Pyx. The derived values of ejecta mass (M e j; top 
panel) and electron temperature (T e ; lower panel) depend on assumed values for the distance to T Pyx and the maximum expansion velocity 
in the nova outburst. Throughout this paper, we assume d = 4.8 kpc and f ma x = 3,000 km s — 1 , as constrained by optical observations; 
these parameters are marked as dashed red lines. 



light curve evolution at all stages. 

By applying observational constraints on two of the 
four fundamental parameters with measurements, we can 
find non-degenerate fits to the other two; for the remain- 
der of our discussion we consider w max and d to be con- 
strained by o ptical spectroscopy a nd light echo imaging 
with the ff5T ([Osborne et al.ll20TTL Sokoloski et al. 2012, 
in preparation). We find that the ejected mass is in the 
range 1.2-3.3 x 10~ 4 M allowing for the 0.5 kpc uncer- 
tainty on the light echo-derived distance and assuming 
the ejecta velocity is in the range 2000-4000 km s — 1 . 
The best fit models for all parameter combinations have 
significant delays in the ejection time relative to the opti- 
cal discovery discovery, showing that this result is driven 
primarily by the early non-detections of the source, and 
is independent of the choice of t> max or d. The thickness 
of the shell, i> m i n /f miK , is largely determined by how flat- 
topped the radio light curve is, and is also independent 
of the assumed values of v max or d. 

4.3. Modeling of Clumpy Ejecta 

In the simple modeling applied thus far, we as- 
sumed that the ejecta are smoothly distributed with 
a volume filling factor of unity. However, there 
are indications from both optical spectroscopy dur- 
ing nova outbursts and imaging of nova remnants 
that ejecta are often clumpy, with fi lling factors in 
the range / = 10~ 5 - 10" 1 (e.g., iWilliamsl U99l 



Andrea et al.lll994HSaizar fc Ferlandlll994l iMason et al.1 
2005t lEderodite et al. 1 120061 : IShara et al.ll2012D . T Pvx 

itself shows significant evidence for clumping in its older 
ejecta, as the surrounding Ha+[N III nebula has been 
resolved into a ring of knots by HST jSh ara et al.1119971 : 
ISchaefer et aT1l27)ll ). 

An observed radio light curve can be produced by 
clumpy ejecta with relatively lower mass and electron 
temperature than is required for homogeneous ejecta. 
We can adapt our simple model of radio emission from 
novae to account for clumping if we assume that a sin- 
gle volume filling factor (/) can describe the ejecta at 
all times and radii, and that that the material between 
clumps is of zero density. We also continue to assume 
that the ejecta are isothermal. In this case, clumping 
simply increases the radio luminosity by a constant fac- 
tor over the uniform-density case, with no changes to 
the time evolution or radio spectr um ([Abbott et al.lll98ll : 
iNugis et al.lll99aiHevwooHl2004D . The radio flux density 
from clumped ejecta is boosted by a factor: 

/1\ 2/3 

""Vclump \~f J *"V uniform- W 

This trend is illustrated in Figure [71 where we hold d and 
w mM fixed and plot how best-fit ejecta mass and electron 
temperature vary with assumed values for /. 
However, if the ejecta in the model are too clumpy, un- 
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Figure 7. Best-fit electron temperatures (upper panel) and 
ejected mass (lower panel) for clumpy ejecta models with a range 
of volume filling factors, assuming d = 4.8 kpc and D max =3,000 
km s" 1 . The red hatched region in the top panel can be excluded 
because the electron temperatures are <5,000 K, i.e. too low for 
hydrogen to be ionized. The vertical dashed line therefore repre- 
sents the minimum filling factor allowable (/ = 0.02). The hori- 
zontal dotted line in the bottom panel marks the ejecta mass that 
corresponds to this minimum filling factor, M cj = 3.7 X 10~ 5 M Q . 



physically low temperatures are required to describe the 
light curve, T e < 5, 000 K. At such low temperatures, hy- 
drogen would not be ionized and no free-free radio flux 
would be emitted. We can therefore exclude the red area 
of parameter space in Figure [7J and place limits on the 
volume filling factor in T Pyx, / > 0.02. This also allows 
us to place a lower limit on the ejecta mass in T Pyx of 
M cj > 3.7 x 10~ 5 M (bottom panel of Figure UJ. Ac- 
counting for the factor of ~2 uncertainty introduced by 
the errors on the distance and ejecta velocity, this lower 
limit can be reduced to ^2 x 10~ 5 M Q . We therefore 
conclude that the ejecta in T Pyx had a total mass in the 
range M C j = (2 — 33) x 10~ 5 M Q , spanning the parameter 
space from maximally clumpy to uniform density. 

5. DISCUSSION 

5.1. Alternative scenarios for producing radio emission 

in T Pyx 

We assert in the previous section that the radio light 
curve is the product of material ejected during the 2011 
nova outburst. In this section we show that other origins 
for the radio light curve can be confidently ruled out. 

Given the extended nebular remnant around T Pyx vis- 
ible in Ha+[N II], it is conceivable that the radio emis- 
sion could originate from this spatially-extended mate- 
rial, as the radiation associated with the outburst inter- 
acts with (and presumably ionizes) the remnant. Indeed, 
the timescales on which we see the radio emission develop 
are appropriat e . Bas ed on the HST imaging presented in 
iSchaefer et al.l (|2010l ). the majority of clumps in the rem- 



nant are located at distances between 2.3 and 5.5 arcsec 
from the central source. Placing T Pyx at 4.8 kpc, these 
angles correspond to physical distances of 1.6-3.9 x 10 17 
cm from the central system. The light crossing times to 
reach the inner and outer extents of the shell are therefore 
62 and 149 days, respectively — corresponding roughly to 
the time we see the radio light curve rise and then fade. 

However, the Ha luminosities observed during the out- 
burst from the extended remnant are much too low for 
the photoionized gas in the nebula to have produced the 
high observed radio flux density. An HST observation 
of T Pyx on 2011 September 25 (Day 164) using the F 
656N filter yields an integrated Ha flux for all knots be- 
yond 3.5 arcsec of ^4 x 10~ 17 erg s _1 cm~ 2 (A. Crotts, 
private communication). Correcting for an interstellar 
E(B-V ) of 0.49 using the extinction law of iFitzpatrickl 
(1999), and using the conversion relation between Ha 
and 5 GHz fluxes expected for a 10 4 K plasma given in 
ICaplan fc Deharven g (1986), we find that the observed 
Ha flux implies a radio flux density at 5 GHz from the 
extended remnant of ~0.1 /iJy — five orders of magni- 
tude lower than observed. Moreover, we never resolve 
the radio source associated with T Pyx, despite the fact 
that the VLA was in its relatively extended A- and B- 
configurations around light curve maximum. A nebula 
with a diameter of a few arcseconds would have been 
thoroughly sampled with >10 synthesized beams across 
it at our higher observing frequencies, but T Pyx was 
consistent with a point-like distribution in all observa- 
tions. 

Another alternative scenario is that the radio emission 
originates in a dense medium in the i mmediate vicin- 
ity of the central binary. For example, Chomiu k et al.l 
(|2012| ) showed that the radio emission in the recent nova 
V407 Cyg originated primarily in the dense wind of the 
Mira giant companion in that system. During the period 
where stable nuclear burning took place on the white 
dwarf's surface, the circumbinary medium was illumi- 
nated and ionized, and the radio flux density from this 
material increased. Could a similar mechanism be at play 
in T Pyx? Long-term optical light light curves of T Pyx 
show a change in the orbital period with time, consistent 
with mass loss rate f rom the companion star as high as 
5 x KT 7 vr" 1 (IPatterson et al.| [l99l lUthas et al.l 
[20101 ISchaefer et alJI2011D . If mass transfer in the sys- 
tem is even mildly non-conservative, then a significant 
amount of material could be lost into the circumbinary 
environment, which might emit in the radio once it is 
ionized by the nova outburst. 

However, radio emission from this circumbinary mate- 
rial should e mit in a manner simila r to a stellar wind (see 
discussion in lChomiuk et aLll2012f) and therefore exhibit 
a partially optically thin spectrum with a « 0.6. The 
spectrum we observe during the rise is more consistent 
with a completely optically-thick body (a « 1.8), and 
the early increase in flux density implies that this body 
is expanding. In addition, if the circumbinary medium is 
in the immediate vicinity of the white dwarf, it would be 
difficult to account for the delay of >50 days between the 
optical rise and the radio rise. Therefore, it is unlikely 
that T Pyx's radio light curve is the result of ionization 
or interaction in a pre-existing circumbinary medium. 

Given these factors, and the relatively good match be- 
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Figure 8. Figure 4 from Nomoto ct al. (2007), plotting accre- 
tion rate against white dwarf mass and showing different regimes 
of accretion. Our estimated range of accretion rates in T Pyx is 
over-plotted as a grey shaded region, assuming that the entirety 
of M e j was accreted since 1966. If instead we eject three times 
more mass than was accreted since 1966, M could be as low as 
the dashed red line. At low accretion rates (below the solid line), 
the accreted material occasionally undergoes thermonuclear run- 
away in a nova outburst; estimated ejecta masses are marked with 
dashed lines. Intermediate accretion rates between the solid and 
dot-dashed result in steady nuclear burning of hydrogen on the 
white dwarf surface. High accretion rates above the dot-dashed 
line lead to unstable configuration wherein the accreted envelope 
expands in a "red giant" -like configuration. 

tween the multi-frequency light curves and the simple 
model of expanding thermal ejecta presented in Section 
[31 we conclude that the expanding nova ejecta are the 
most likely origin of the radio emission. 

5.2. Implications of a massive ejection event in T Pyx 

The ejected masses derived from our model fits to the 
radio light curve are large, ranging from 2 x 10 -5 M Q 
for clumped ejecta, to 3.3 x 10 -4 M for a shell with a 
smooth density distribution of the for n oc r~ 2 . Despite 
some uncertainties in the geometry, density profile, and 
temperature evolution of the ejecta, the primary con- 
straint on the ejecta mass is simply provided by the high 
flux density observed at light curve maximum. We note 
that lSelvelli et all (|2008l ) estimated an ejecta mass in the 
range 10 -5 -10 -4 for the 1966 outburst based on the long 
duration of the optically thick phase in the optical light 
curve — the radio data presented here support that con- 
clusion. This high ejected mass challenges commonly 
held assumptions about T Pyx and holds important im- 
plications for the global parameters of the system. 

Since we know the recurrence time of T Pyx, we can es- 
timate the accretion rate in the system. Assuming that 
all the material ejected in the 2011 outburst has been 
accreted since the 1966 outburst, the average accretion 
rate is simply given by M = For a recurrence time 
tree = 44 years, our estimated range of M e j implies an 



accretion rate in the range M — 4.5 x 1CP 7 — 7.5 x 
10 -6 M Q yr" 1 . Theoretical studies of nova outbursts 
find that some novae eject as m uch as three times the 
accreted mass (jYaron et al.l 120051 ) . in which case the ac- 
cretion rate could be up to a factor of ~3 lower (i.e. 
1.5 x 10~ 7 - 2.5 x 10~ 6 M Q yr- 1 .) 

Estimates of the accretion rate from observations im- 
ply accumulated masses since the 1966 outburst that are 
only just compatible with the very lowest ejecta masses 
we have derived from the radio data — and only if the 
nova outburst ejects more material than wa s previously 
accre ted. Assuming E(B-V) = 0.49 ± 0.17 ([Shore et all 
120111) . and a distance to T Pyx of 4.8 kpc, the UV flux 
reported in lSelvelli et al.l (|2008D is consistent with a mass 
accret ion rate as high as ~ 10 _7 M Q yr -1 . Uth as et all 
(2010) present a new determination of the orbital period 
of T Pyx based on 14 years of photometric data, and find 
highly significant evidence of a period slow down (con- 
firming the earlier result of iPatterson et ail (|1998D ) that 
implies mass loss from the companion star. They derive 
a mass loss rate of 5 x 10~ 7 M Q yr -1 using their estimate 
of the mass ratio in the system, and note that this rep- 
resents the upper limit to the mass accretion onto the 
white dwarf — the rate could be lower if mass transfer is 
non-conservative. 

The range of accretion rates implied by our radio light 
curves is also uncomfortably high in light of theoretical 
expectations for the fate of accreted material on a white 
dwarf. Acco rding to studies of thermal stability on white 
dwarfs (e.g. JNomot^ll982t iNomoto et alJ 120071 ). at low 
accretion rates (M < 10~ 7 M Q yr -1 ), accreted material 
collects on the white dwarf until a critical pressure P cl -it 
is reached and a nova outburst is triggered (Figured]). In 
a relatively narrow range of higher accretion rates {M w 
[1 — 6] x 10~ 7 Mq yr -1 ), the white dwarf can stably fuse 
hydrogen in the accreted shell. Finally, at the highest 
accretion rates (M > 6 x 10~ 7 M© yr -1 ), it is likely 
that mass transfer onto the white dwarf will be unstable, 
and the accreted material will puff up into an extended 
red-giant-like atmosphere, or be expelled from the white 
dwar f in an optically-thick accretion wind ([Hachisu et all 
ri999h . 

T Pyx produces novae despite experiencing an accre- 
tion rate that is ~l-2 orders of magnitude higher than 
predicted for nova-hosting white dwarfs. The accretion 
rates implied by our radio-derived ejecta masses and the 
recurrence time in T Pyx are so high that T Pyx was ex- 
pected to have been in an unstable red-giant-like configu- 
ration during quiescence (grey bar in Figure |5J). Even un- 
der the most conservative assumptions of the clumpiest, 
lowest-mass ejection allowed by our radio light curves, 
and an ejected mass that is three times larger than the 
accreted mass, the accretion rate is still high enough that 
the white dwarf would have at least been stably burn- 
ing hydrogen on its surface in the standard theoretical 
picture. If the white dwarf had been in such a sta- 
ble, hydrogen-burning configuration during quiescence, 
it would not have undergone a nova outburst in 2011. 
Indeed, a hypothes is that T Pyx might host a super- 
soft X-ray source ([Patterson et alJ [l998: Knigg e et all 
2000) was over-turned by X-ray observations during qui- 
escence, which showed no evidence for stable nuclear 
burning (|Selvelli et al.l l2008) . In addition, there is no ev- 
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idence for a red-giant-like exten ded envelope or optically- 
thick wind during quiescence (|Girmozzi k Selvellil 120071 ; 
ISelvelli et al.l[2008h . T Pyx may therefore present a sig- 
nificant challenge to the theory of accretion onto white 
dwarfs, and war rant a re-investigatio n of the parameter 
space laid out bv lNomoto et al.1 ()2007l ) . as plotted in Fig- 
ure [51 Our estimate for the average rate of mass transfer 
in T Pyx can only be brought into line with accretion 
theory if we allow the distance to be significantly lower 
than the 4.8 ± 0.5 kpc assumed here, and let the filling 
factor and T e of the ejecta to take on their minimum 
allowable values. 

The large ejecta mass measured during the 2011 out- 
burst of T Pyx implies that the white dwarf has a mass 
that is lower than usual for a recurrent nova. Thermonu- 
clear runaways on the surface of accreting white dwarfs 
are triggered when the pressure in the accreted shell 



GM WD m acc 



(2) 



reaches P, 



w 10 19 dyn cm" 2 (|Fuiimotd 119821 : 
lMacDonaldlll983l ). Assuming that the ejected mass and 
the mass accreted between nova outbursts are essentially 
equivalent, and using the relationship between whit e 
dwarf mass and radius from|Hamada k Salpeterl (|1961l ). 
we find that white dwarfs with masses between 0.65 and 
1.1 M reach P cr it with our estimated range of M e j. If 
the nova ejected three times more mass than it accreted, 
the white dwarf mass could be as high as 1.2 M©. For 
higher mass white dwarfs, a nova outburst would have 
been triggered at a significantly lower accreted mass. 

The results of detailed modeling of accretion-outburst 
cycles in novae confirm this b asic reasonin g . In t he full 
grid of models presented by lYaron et al.l ((2005), only 
white dwarfs with masses <1.25 M ffi can ever ac c rete a 
shell of 10" 5 M Q or greater. iTownslev k Bildstenl (|2004fl 
note that high mass accretion rates, such as the one we 
have inferred in T Pyx, actually reduce the TNR trigger 
mass for a given white dwarf mass relative to low accre- 
tion rates by up to an order of magnitude. The mass esti- 
mates we have calculated using Equation 2 are therefore 
only u pper limits to the white dwarf mass. lUthas et al.l 
(2010) find a white dwarf mass of 0.7 ± 0.2 M based 
on time-resolved optical spectroscopy and assumptions 
about the density of the secondary. Although the deriva- 
tion of white dwarf mass from a spectroscopic period 
depends on several uncertain assumptions, the data are 
compatible with a low-mass white dwarf. Higher white 
dwarf masses (generally 1.3 M ) are quoted extensively 
in the literature for T Pyx, although this value appears 
to be based primarily on the expectation of a high-mass 
white dwarf in system wit h short nova recurrence time 
(e.g.. iWebbink et aLlll987l ). Our results present a signif- 
icant challenge to this assumption. 

5.3. A prolonged, multi-stage ejection in T Pyx 

The radio evolution we have observed in T Pyx is not 
consistent with a single impulsive ejection of the bulk of 
material on day zero. In fact, there is eviden ce that mass 
ejectio n took place in at least two stages. iShore et"aLl 
(|2011[ ) derive an ejected mass from optical spectroscopy 
obtained on Days 1 and 9 of M e j «2x 10~ 5 / M Q , where 
/ is the volume filling factor. This is incompatible with 



our radio derived mass range, since inconsistent filling 
factors are required to bring the estimates into agree- 
ment. However, it is possible that these spectra trace 
only the earliest epoch of mass ejection. We note that 
the single 33 GHz detection of T Pyx with the VLA on 
Day 14 is consistent with optically thick emission, and 
therefore could originate in an ejected shell. The low 
flux density, and subsequent radio non-detections indi- 
cate that the mass associated with this ejection must 
small, at mos t 10~ 6 M m . This is consistent with the 
estimate from IShore et al.l (|2011l ) if the filling factor is 
small, / <0.1. 

The bulk of the ejection of material in T Pyx did not 
happen until several months after optical discovery. As 
we outlined in Section 4.1, the late rise of the radio emis- 
sion indicates that the primary mass ejection took place 
50-100 Days after to. However, even allowing for a late 
ejection, the light curves deviate from the S v oc t 2 rise 
expected from an optically thick shell expanding at con- 
stant velocity. The observed flux density increase be- 
tween Days 66 and 194 cannot be reconciled with ejecta 
expanding at a constant velocity emitted around Day 66. 
In fact, the C-band light curve observed between Days 
45 and 102 rise more as S v cx t , which may indicate a 
more gradual ejection of material, although the data are 
insufficient to draw a definitive conclusion. 

Regardless of the exact details of the earliest stages of 
mass ejection traced by the main radio rise, it is clear 
that T Pyx experienced a significant delay in the ejec- 
tion of the bulk of its accreted shell. This delay may pro- 
vide an explanation for the optical evolution of T Pyx. 
As evidenced in the AAVSO optical light curve shown 
in Figure 1, T Pyx shows a dramatic plateau in optical 
brightness, and does not begin its optical decline phase 
until ^90 days after the start of the outburst. The time 
of this transition in the optical light curve roughly coin- 
cides with the rise of the radio light curve, as expected 
for an optical photosphere that has been approximately 
constant in size for a few months, and suddenly begins 
to recede through the ejecta as they start expanding. 

Other novae have shown similar plateaus around op- 
tical maximum, with the duration of the plateau widely 
varyi ng between sources, fr om a few to hundreds of days 
(e.g. , iKato k H achisu 2 0111) . The physical origin of these 
plateaus around maximum is uncl ear, although the y gen- 
erally last longer on slower novae. iFriedjungl ()1992l ) noted 
that during the plateau in HR Del, the optical spectrum 
showed no evidence of an optically-thick wind as typi- 
cally seen in other novae; it was in fact was more consis- 
tent with an almost stationary atmosphere. The author 
proposed that the nova outburst in this system was ini- 
tially not strong enough to eject the accreted envelope, 
and instead settled into a quasi-static configuration with 
little mass loss. No explanation was offered, however, 
for the transition to the subsequent phase of significant 
mass loss. The plateau around maximum in T Pyx could 
be interpreted in the same way, and indeed the late ejec- 
tion implied by the radio light curves lends credence to 
the idea of a stable, quasi-hydrostatic phase in the early 
stages of outburst in T Pyx. If the nova had been losing 
significant amounts of mass at high velocities from the 
onset of outburst, we would have detected the source at 
h igher radio flux d en sities during our early observations. 

lHachisu k Katol (|2004l ) concurred that plateaus 
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around optical maximum could be interpreted as a static 
phase of the outburst, with an inflated but stable white 
dwarf photosphere formed after the o nset of the ther- 
monu clear runaway. In later work (jKato &; Hachisul 
l20ll . the same authors proposed that the friction pro- 
vided by the companion star interacting with the hydro- 
static envelope could provide an additional source of en- 
ergy, empowering the outburst to transition from a stable 
to a mass-losing configuration. Although we do not spec- 
ulate as to whether orbital friction played a role in the 
observed late ejection in T Pyx, it is worth noting that 
iKato fc Hachisul ([20111 ) find that this mechanism should 
be most efficient in very short period binaries, like T Pyx, 
where the companion becomes deeply embedded in the 
nova envelope. 

6. CONCLUSIONS 

From our analysis of multi-frequency radio light curves 
obtained with the newly-upgraded VLA, we reach several 
conclusions that challenge some of the most fundamen- 
tal and commonly-held beliefs about the T Pyx binary 
system. 

1. The 2011 nova event in T Pyx expelled a large 
amount of mass, totaling [2 — 33] x 10~ 5 M©. 

2. In order to have accreted sufficient material since 
1966 to account for this large ejecta mass, the white 
dwarf in T Pyx must have been accreting material 
at a very high rate, M > 1.5 x 10~ 7 Af© yr" 1 . 

3. Such a high nova trigger mass implies that the 
white dwarf in T Pyx is not near the Chan- 
drasekhar mass. Instead, it more likely has a mass 
of <1 M . 

4. The bulk of the mass was not ejected from the 
white dwarf until ~50-100 days after the optical 
discovery. 
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Table 1 

VLA Observations of T Pyx 



Date 
CUT) 


t-t a 
(Days) 


Epoch 


Conf. 


Freq 
(GHz) 


S v 
(mjy) 


Freq 
(GHz) 


S v 
(mjy) 


Freq 
(GHz) 


(mjy) 


Freq 
(GHz) 


S u 
(mjy) 


2011 Apr 22.0 
2011 Apr 22.2 


8.0 
8.2 


1 
1 


B 
13 


5.9 
33.1 


0.003 ± 0.008 
0.016 ± 0.030 














2011 Apr 30.9 
2011 May 1.1 


16.9 
17.1 


2 
2 


13 
13 


5.9 
33.1 


0.027 ±0.011 
0.248 ± 0.040 














2011 May 15.5 
2011 May 15.9 


31.5 
31.9 


3 
3 


BnA 
BnA 


1.6 
3.0 


-0.040 ± 0.031 
0.000 ±0.035 


6.0 

8.4 


0.022 ± 0.006 
0.017 ±0.018 


14.3 
22.0 


0.028 ± 0.017 
0.038 ± 0.049 


33.0 


-0.046 ± 0.055 


2011 May 29.0 


45.0 


4 


BnA 


1.5 
14.3 


0.134 ±0.068 
0.026 ± 0.019 


3.0 
22.0 


0.015 ±0.034 
0.021 ± 0.031 


6.0 
33.0 


0.005 ±0.010 
0.011 ± 0.030 


8.4 


0.032 ±0.017 


2011 Jun 1.0 


48.0 


5 


BnA— >A 


33.1 


0.048 ± 0.029 














2011 Jun 19.9 


66.9 


6 


A 


5.0 


0.130 ±0.017 


7.0 


0.194 ±0.016 










2011 Jul 24.8 


101.8 


7 


A 


1.6 
7.0 


0.082 ± 0.040 
0.529 ±0.039 


2.5 


0.119 ±0.052 


3.5 


0.185 ±0.039 


5.0 


0.304 ± 0.031 


2011 Aug 20.7 


128.7 


8 


A 


1.3 
5.0 
16.0 


0.146 ±0.058 
1.128 ± 0.061 
6.932 ±0.715 


1.8 
7.0 

22.0 


0.190 ±0.066 

I. 885 ± 0.097 

II. 38 ± 1.17 


2.5 
8.4 
29.0 


0.318 ±0.056 
2.809 ± 0.141 
18.82 ± 1.91 


3.5 
12.6 
37.0 


0.581 ±0.065 
4.832 ± 0.500 
23.28 ± 2.36 


2011 Sep 17.6 


156.6 


9 


A^D 


12.6 
37.0 


11.09 ±0.12 
47.90 ± 4.94 


16.0 


14.20 ± 1.44 


23.7 


26.93 ± 2.71 


29.0 


38.97 ±3.94 


2011 Sep 24.6 


163.6 


10 


A^D 


2.5 

8.4 


0.800 ±0.126 
6.738 ±0.338 


3.5 


1.350 ±0.098 


5.0 


3.014 ±0.161 


7.0 


5.157 ±0.260 


2011 Oct 25.5 


194.5 


11 


D 


1.3 

5.0 
23.7 


-0.763 ±0.210 
3.816 ±0.196 
26.27 ±2.63 


1.8 

7.0 
29.0 


0.631 ±0.210 
6.106 ±0.308 
30.97 ±3.10 


2.5 
12.9 
37.0 


1.421 ±0.139 
14.09 ± 1.41 
34.68 ± 3.47 


3.5 
16.0 


1.716 ±0.127 
18.60 ± 1.86 


2011 Nov 10.5 


210.5 


12 


D 


1.3 
5.0 
23.7 


0.057 ± 0.243 
3.685 ±0.190 
26.32 ±2.63 


1.8 

7.0 
29.0 


0.774 ± 0.243 
6.263 ±0.315 
30.51 ± 3.05 


2.5 
12.6 
37.0 


1.166 ±0.136 
14.16 ± 1.42 
33.40 ± 3.35 


3.5 
16.0 


1.797 ±0.124 
17.64 ± 1.77 


2011 Dec 7.4 


237.4 


13 


D 


1.3 
5.0 
23.7 


0.699 ±0.277 
4.270 ±0.219 
18.14 ±1.82 


1.8 

7.0 
29.0 


0.445 ±0.167 
6.594 ±0.332 
18.87 ± 1.89 


2.5 
13.3 
37.0 


1.567 ±0.121 
12.41 ± 1.25 
19.80 ± 1.99 


3.5 
16.0 


2.307 ±0.128 
14.37 ± 1.44 


2012 Jan 5.4 


266.4 


11 


DnC 


1.3 
5.0 
23.7 


0.277 ±0.210 
4.345 ± 0.222 
10.12 ± 1.01 


1.8 

7.0 
29.0 


0.477 ±0.121 
6.096 ± 0.307 
9.449 ±0.952 


2.5 

13.3 
37.0 


1.528 ±0.146 
8.995 ± 0.903 
9.061 ±0.926 


3.5 
16.0 


2.659 ±0.153 
9.434 ± 0.947 


2012 Feb 8.2 


300.2 


15 


C 


1.3 
5.0 
23.7 


1.000 ± 0.153 
3.802 ±0.193 
5.155 ±0.519 


1.8 

7.0 
29.0 


1.441 ± 0.161 
4.685 ±0.236 
5.045 ±0.511 


2.5 
13.3 
37.0 


1.876 ± 0.137 
5.596 ±0.561 
4.385 ±0.454 


3.5 
16.0 


2.580 ± 0.143 
5.739 ±0.575 


2012 Mar 10.1 


331.1 


16 


c 


2.5 
13.3 
37.0 


1.757 ±0.112 
3.981 ± 0.401 
2.826 ±0.329 


3.5 
16.0 


2.335 ±0.127 
4.125 ± 0.416 


5.0 
23.7 


3.085 ±0.160 
3.645 ± 0.371 


7.0 
29.0 


3.492 ±0.177 
3.031 ± 0.324 


2012 Mar 27.1 


348.1 


17 


c 


1.3 


0.529 ±0.176 


1.8 


0.927 ±0.125 










2012 Jul 22.7 


465.7 


18 


B 


1.4 
5.0 


0.600 ±0.133 
0.867 ±0.075 


1.8 
7.0 


0.592 ±0.098 
0.812 ±0.074 


2.6 


0.764 ±0.104 


3.5 


0.689 ± 0.070 


2012 Sep 07.6 


512.6 


19 


BnA 


1.4 
5.0 


0.622 ±0.132 
0.669 ± 0.065 


1.7 
7.0 


0.562 ±0.077 
0.496 ± 0.048 


2.6 


0.683 ±0.083 


3.5 


0.642 ±0.058 


2012 Sep 23.5 


528.5 


20 


BnA 


13.3 


0.757 ±0.096 


17.5 


0.682 ±0.104 


23.7 


0.866 ±0.214 


29.0 


0.483 ±0.151 



37.0 0.464 ±0.182 



a Time since optical discovery of the nova event (2011 April 14; Waagan ct al. 2011). 



